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I. 


Background 

As  material  for  possible  historical  (as  well  as  motivational) 
interest,  we  may  observe  that  Volume  1,  No.  1  of  Management  Science 
contained  an  article  entitled  "The  Stepping  Stone  Method  for  Explaining 
Linear  Programming  Calculations  in  Transportation  Type  Models.  " 
Volume  1,  No.  2  contained  an  article  entitled  "Optimal  Estimation  of 
Executive  Compensation  by  Linear  Programming.  "  The  purpose  of 
these  two  papers  was  to  help  this  new  journal  and  the  new  society, 

TIMS,  which  published  it  get  off  to  a  start  that  would  also  provide  a 
basis  for  further  growth  not  only  in  themselves  but  also  in  a  society 
(TIMS)  which  would,  in  turn,  greatly  enlarge  the  prospects  for  these 
and  other  related  developments  in  the  management  sciences. 

By  emphasis  at  least,  the  first  of  these  two  articles  was  designed 

to  appeal  to  immediate  use  and  understanding.  The  emphasis  in  the 

second  article  was  directed  to  longer  range  potentials  for  additional 

applications  in  the  then  new  discipline  of  linear  programming.  Some 

of  these  applications  included  extensions  to  areas  like  "inequality  con- 

strared  regression,"  "htfulti-1  dimensional  objective  optimizations" 

and  their  extensions  to  ordinal  non-metric  scaling,  including  non- 

2 

Archimedean  constructs.  New  theory  as  well  as  new  methodological 
possibility  s  were  also  haturally  kept  in  mind,  and  exploited  as  their 


1.  See  1 2]  and  [3]  in  the  bibliography  appendix  to  this  paper. 

2.  See  (4]. 
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potenoy  for  applications  appeared  to  warrant.  Thus,  the  transportation 
model  in  the  first  of  these  two  articles  was  extended  via  a  variety  of  new 

n 

ideas  which  ranged  from  the  development  of  the  "poly  -  m  method'0  to  the 

use  of  model  transformations  and  approximations.  Ideas  in  the  second 

article  subsequently  gave  birth  to  attempts  to  exploit  the  use  of  linear 

programming  models  and  methods  to  attack  a  variety  of  nonlinear 
4 

problems.  They  also  gave  rise  to  the  gO®l  programming  ideas  which 
have  now  begun  to  be  used  in  a  variety  of  novel  contexts.  One  such  use, 
the  OCMM^series  of  models  [121,  has  also  involved  a  change  in  the 
state  of  modeling  per  se  by  joining  the  ideas  of  goal  programming  and 
Markov  processes  into  a  new  approach  for  dealing  with  decision  (Markov) 
processes. 

As  in  these  other  developments,  applications  associated  with 
the  OCMM  models  have  provided  points  which  now  make  it 
possible  to  undertake  further  methodological  and  theoretical  developments, 
as  well  as  applications,  that  will  provide  still  better  bases  for  future 
extensions.  This,  in  fact,  is  the  point  of  the  present  paper,  which  will 
be  followed  by  others,  in  which  we  shall  undertake  to  join  some  of  the 
preceding  work  (e.  g. ,  in  model  approximation  and  goal  programming) 

3.  See  [4]. 

4.  See  fill,  1131.  [14]. 

5.  OCMM  (Office  of  Civilian  Manpower  Management,  U.  S.  Navy) 
This  series  of  models  is  described  in  [12]. 
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not  only  with  each  other  but  also  with  presently  disparate  developments 
(e.g. ,  in  linear  programming  under  uncertainty)  while  we  simultaneously 
attempt  to  provide  solution  procedures  (including  good  starts)  that  can 
supply  operational  significance  for  the  new  problem  areas  that  might 
then  be  addressed. 

Introduction 

Explicit  solutions  to  linear  programming  or  linear  fractional 

programming  problems  presently  exist  only  for  special  cases,  the 

6 

most  general  of  which  are  the  full  row  rank  interval  programming 

7 

and  the  full  row  rank  linear  fractional  programming. 

In  some  cases  these  models  are  of  direct  interest.  This 

naturally  includes  those  cases  where  problems  can  be  modelled  to 

fit  the  theoretical  and  computational  requirements.  It  also  includes 

cases  for  which  model  (structurd  approximation  and  parameterization 

methods  can  yield  either  exact  or  approximate  solutions  to  more 

8 

complicated  problems.  Finally  it  includes  cases  in  which  advanced 
starts  can  be  thus  obtained  including  ones  which  greatly  accelerate 
the  attainment  of  solutions--as  when,  e.g.  ,  primal  simplex  technique 
applied  to  such  starts  have  yielded  great  improvements  in  efficiency 

6.  See  |1). 

7.  See  [  5j. 

8.  See  [  6l. 
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relative  to  alternative  solution  procedures  that  might  otherwise  have 
been  required. ^ 

There  are  at  present  no  comparable  classes  of  explicit 

solutions  to  "goal  programming"  models.  The  increased  usage  now 

being  accorded  to  this  type  of  model,  however,  would  seem  to  justify 

such  developments  in  order  to  anticipate  benefits  for  this  class  which 

are  analogous  to  those  described  in  the  preceding  paragraph. 

The  present  paper  is  directed  toward  developing  such  solution 

procedures  for  the  general  case  of  piecewise  linear,  but  separable, 

"goal  functionals.  "  We  define  a  separable  functional  as  a  sum  dt 

one -variable  goal  functionals.  A  one -variable  goal  functional,  f.(x  ), 

J  J 

is  a  function  that  is  monotone  decreasing  for  x  *  g  ~-g.  *  the  stip- 

J  j  3 

ulated  jth  goal- -and  monotone  increasing  for  x-  >g..  Note,  however, 

J  J 

that  strict  monotonicity  is  not  required.  In  fact  for  the  important 
"goal  interval"  type  of  model,  the  functional  has  a  constant  minimum 
over  a  whole  interval. 

The  following  examples  are  included  in  this  class: 

(i)  Absolute  value  functions,  including  those  with  asymmetric 
weights  and  multi-goal  components. 

(Ci )  General  convex  piecewise  linear  functionals 

(iii)  Goal  interval  functionals 


9.  See  (8)  and  [10). 
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(i v >  Ilypermodian  functionals  and  related  functionals  in  extensioi 
of  ordinary  goal  programming.  ^ 

This  list  is  not  exhaustive,  but  it  is  nevertheless  indicative 
of  the  manv  significant  kinds  of  problems  that  are  included  for  the 
development  that  we  shall  now  provide  as  follows:  First  we  shall  state 
and  prove  a  very  general  lemma  which  can  be  used  to  reduce  and 
simplify  such  problems.  This  lemma  will  be  presented  in  the  next 
section--i.  e.  ,  section  2--along  with  its  application  to  a  simple  (but 
significant)  class  of  cases.  The  lemma,  we  may  note,  is  not  restricted 
to  the  latter  class  of  cases  and  is,  indeed,  established  in  far  greater 
generality  than  required  for  any  part  of  the  present  paper.  In  addition, 
we  shall  proceed  to  our  results  for  obtaining  explicit  solutions  via 
general  separable  convex  piecewise  linear  functionals,  which  will  be 
covered  in  section  3.  In  this  section--i. e. ,  section  3--we  shall 
obtain  a  linear  programming  equivalent  which  will  be  extended  in 
section  4  where  procedures  for  obtaining  explicit  solutions  will  be 
delineated.  Section  5  will  then  conclude  tti  s  paper  with  a  numerical 
example  which  will  serve  to  illustrate  some  of  the  results  achieved 
in  the  present  paper  and  also  indicate  some  of  the  possible  extensions 
that  will  be  covered  in  later  papers  in  this  series. 


10.  Sec  [  7  ] 
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2.  A  General  Reduction  Lemma: 


The  following  lemma  is  established  in  far  greater  generality 

than  we  shall  require  for  this  paper.  Indeed,  no  properties  other 

than  monotone  decrease  up  to  the  goal  and  monotone  increase  to  its 

right  are  stipulated.  Thus,  we  let  f  =Ef  (x.)  where  f.(x  )  is  monotone 

j  j  J  J  J 

decreasing  for  Xj  sh^  and  monotone  increasing  for  2  h^.  Then  we 
consider  the  following  problem: 

(1.1)  Min.  f(x )  =  £f.  (x  ) 


J 


J  J 


subject  to  £x,  sb. 

3  J  1 

(1.  2)  a_,  s  Xj  «b. 

With  this  in  mind  we  now  develop 


Lemma  1:  If  x  satisfies  (1.2)  and  for  some  j0,  x-  >  h,  ,  then 
5  satisfies  (1.2)  and  f(y)  s  f(x),  where  °  ° 

x  ■  (.  . . ,  JTj,  . . . ), 


x  =  f  Xj*  ^  }t> 

3  \jnax  (ajo,  hj(j),  j=j0. 


Proof:  If  x^  is  replaced  by  xJq  -  A  ,  Os  As  min.  (xJq-  h^  ,  Xjo-aJo),  then 

°  L  (x.  -A)sr  (x.  ) 

Jo  3o  J0  Jo 

a  s  X;  -  As  b. 

J  Jo  3 

and  £x.  +  x  -A  s  b  ,  where  j  signifies  that  the  summation 


Jo 


3  Jo 


omits  the  variable  with  this  index. 
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Thus,  (1.2)  is  still  satisfied.  Also  f  (x,  )  is  monotone  increasing 

Jo  Jo 

for  h;  *  x.  .  We  therefore  obtain  a  decrease  (or  at  lease  no 
Jo  Jo  _ 

increase)  by  setting  A  at  its  upper  limit  S’  and  define  XjQ  =  Xj  -  A 


max  (h*  ,  a4  )  s  x,  . 
Jo  Jo  Jo 


But  x4  >  h4  .  Hence 
Jo  Jo 


and 


W  s  fJo<xJo’ 
£fj(xj,+  ri0WIo,S,(xl- 


I.  e.  .  J 

f(x)  sf(x), 

and  the  lemma  is  thus  established. 


The  lemma  covers  in  particular  nonlinear  goal  functionals  of  types 
such  as  those  drawn  in  the  following  figures. 


Non-Convex  Functional 
Figure  la 


Convex  Functional 
Figure  lb 


Piecewise  Linear 
Functional 
Figure  lc 


Note,  further,  also  that  the  simple  constraint £)Xj  sbQ  in  the  lemma  may 
be  replaced  by  any  system  of  constraints  such  that  decreasing  any  single 
variable  preserves  feasibility  of  the  constraints.  E.  g.  ,  the  system 
5xj  +  2x„  +  x4s  20 

x2  +  3xM  +  5x4  s50 
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would  be  a  valid  type  of  replacement  for  the  single  constraint  of  the  lemma. 
See  the  illustration  in  Section  5. 

We  proceed  first  to  consider,  however,  the  particular  case 

(2. 1)  min  £  U.  |x,  -  g  | 

J  1 

subject  to  £x .  s  bQ 

j 

(2.21 

where,  without  loss  of  generality,  the  constants  Uj  are  indexed  so  that 

Mj  *  2  . . .  *  .  .  .  2  0.  Because  of  Lemma  1,  we  can  impose  the 

additional  condition 

(3.0)  x.  *  max  (g  ,  a^sb'.. 

J  J  J  J 

Via  (2.  2),  however,  a.sx.  so  that  if  g4  *a.  we  need  only  set  x.=a,  and 

3  J  J  J  1  1 

reduce  the  preceding  problem  to  one  in  the  remaining  variables 

(3.1)  minE  U.  |  x.  -  g.| 

j  j  •>  j 

subject  to 

2Vb'° 

J 

(3.  2)  ai  *  x,  s  b  i 

J  J  J 

*  -  / 
where  b  ^  is  defined  as  in  (3.  0)  and  a^  *  g^  £  b  j  with,  again,  Mj  2  U2  •  • * 

Mj  2  a  Mn  i  0,  by  re  numbering,  if  required,  and  b^  is  bQ  reduced 

by  the  x  =  a=  values.  Abusing  notation,  we  replace  n  ”  by  n  too. 

J  J  , 


Because  of  Lemma  1  and  the  definition  of  b  ,  however,  we  now 

j 

have  xj  "  gj  all  3  and  this  implies 


lxj-g3,  =  grxj 
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Thus  w(  may  now  replant  the  preceding  problem  by 


(4.1)  max x 

.11  I 

subject  to  i->  x.  sb' 


j  3 


(4.2) 


a. 


3  J 

This  is  an  ordinary  linear  programming  problem  with  a  dual  that  mav  be 


written 


(5.1)  minvobo  +  E(v  b.'-  v  aJ 
u  °  j  J  J  J  J 

Subject  to 


(5.2) 


+  v  -  v 
J  J 

V  V  vi  *  °- 


=  H 


If2ja  >5^  then  the  constraints  of  (4.  2)  are  evidently  inconsistent 
J  J 

since  the  a  are  lower  bounds  which  the  corresponding  x  must  satisfy. 

•I  J 

If  we  assume  that£a:*b0  then  the  optimum  solution  of  (4.  2)  can  be 

j  J 

immediately  written: 

xi  =  j*l . k-1 

x’  -  t/  -  -  £  a 

k  o  j =1  J  j=k+l  J 

x.l  =  aj*  J  =  k  +  l'  ’  n 

where  k  is  the  smallest  positive  integer  such  tha*. 
k-1  ,  ^  n 

O.  b:  £  bn  ”  £  S: 

.1  =1  3  °  k+1  J 


b >  b *  ~ VT  3 
J*1  J  0  K^l  3 


and 
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To  see  that  these  x*  choices  are  optimal  we  observe  thayhey  satisfy 

F  W  x*  b^  +  W  (b*-£b  -  Zj  a.  + 

the  constraints  in  (4,2)  and  produce  j  j  1  J  J  kojJ  k+1  J  ^-t)  i  j 
in  the  functional  (4.  1).  Via  regrouping11  this  becomes  +  !^  (Mj_wk)bjk*l<Ml "  Mk,a 

Applying  the  "regrouping  principle"  (see  [  )),  we  posit  these  coefficients 
as  optimal  dual  variables  and  need  only  verify  that  they  satisfy  the  dual 
constraints  -  -viz., 


=  M, 
o  k 


k-1 


V  V^*°’  j-1. 

v*-  =  o,  j=l . k-1 


vk+  =  v*"  =  ° 

v*"»  Mk  '  a0>  ,  n 


v  =  0,  j=  k-H. 


n. 


■H 


These  values  of  vQ,  Vj,  v^  are  non-negative  as  required.  They  also 
satisfy  the  other  dual  constraints  in  (5.  2)  since,  as  may  be  seen,  these 
v*  choices  give 

Mk  +  Uj-Uk^.  i=l . k-1 

Mu  =  Mu 


Wk  +  ui  '  uk  s  wj* 


Jrk+1, . .  .  ,  n 


upon  substitution  in  (5.2). 


11.  I.  e. ,  we  are  here  using  a  variant  of  the  "regrouping 
principle"  first  set  forth  in  f  9  ]  as  a  means  of  obtaining  explicit  solutions 
via  the  dual  problem.  See  also  (4|. 


ed3 
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This  completes  the  proof  that  our  explicit  solution  x*  is  optimal. 
It  also  shows  in  full  detail  how  we  can  use  the  lemma  to  obtain  explicit 
solutions  in  this  class  of  cases.  The  next  section  will  extend  this 
result  to  the  case  of  general  separable  convex  piecewise  linear 
functionals. 

3.  The  General  Separable  Convex  Piecewise  Linear  Functional 

The  functional  we  are  considering  is  separable--i.  e. ,  it  is 
a  sun.  of  functions  each  involving  only  one  variable.  Hence,  we  will 
be  able  to  apply  Lemma  1  and  perform  a  reduction  on  one  variable 
at  a  time  to  obtain  a  linear  programming  format.  The  explicit 
solution  will  then  be  developed  and  exhibited  in  the  section  after  this 
one. 

The  general  (continuous)  one  variable  case  has  a  graph  con¬ 
sisting  of  straight  line  segments  with  slopes  that  increase  from  left 
to  right.  See  Figure  lc.  As  is  well  known  classically,  such  a  function 
can  be  written  as 

(6.1)  fa(t)  |  t-g.  |+  pt  +  q 

where  the  g  ,  M  ,  p  and  q  are  given  constants  with  M  =  0,  all  j. 

.1  J  •* 

In  this  form  we  state  the  following: 

Lemma  2:  A  necessary  and  sufficient  condition  for  fa(t)  with 
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Another  representation  of  such  a  function  is  also  possible,  bearing  in  mind 

that  we  shall  wish  to  minimize  on  fit).  Because  it  explicitly  specifies  the 

slopes  involved  in  fit)  we  shall  specify  It  by  f8(t).  We  shall  now  proceed 

to  develop  this  second  representation  of  fit)  and  then  relate  it  to  the  f^it) 

representation  that  we  have  already  given. 

Without  loss  of  generality,  we  may  suppose  that  g^  s  g^  +j  for  all 

j.  Suppose  also  that  a  *  t  <  b.  We  may  then  write  t  =  2t.,  for  some 

J  • 

sequence  of  F  starting  with  t^,  wherein  we  restrict  the  F  by  the  following 
conditions : 


*1*81 

0  ‘  *  gj  '  *j-l*  j*2 . n 

gns  Vl 

Since  we  are  minimizing  on  fit)  for  a  f  t  «  b,  we  will  choose  t  >  0  only  if 

J 

tj-l  (which  is  associated  with  the  greater  slope)  has  reached  its  upper 

bound.  Thus,  for  minimization  purposes,  f(t)  convex  can  be  equivalently 

written  in  separated  linear  form  as 

(6.2)  fs(t)  =  2,  M  +  d, 

J=1  1  J  1 

where  kj  is  the  slope  of  the  line  segment  of  the  graph  of  f(t)  which  runs 
fron>  g(_j  to  g(. 

For  various  important  goal  functionals  we  shall  need  to  be  able  to 
pass  from  one  representation  of  the  function  to  the  other.  We  proceed  to 
develop  the  needed  formulae. 
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Consider  then 


and 


ftt(t)  =£  (i.  I  t-g  |  +  pt  +  q 
j-l  J  J 


_  n+1 

fs(t)  =  £k  U  td,  k:  sk..,,  Vj. 


j=l 


j  J  *  W 


The  graph  of  fH)  is  shown  in  Figure  2.  Now  consider 

*1**1 

Oct  *g  -  g  i 

j  V* 


Oc  t  . 
n+1 

Fort+®:  fa(t)  +  P>Vl 

f8(t)  -  Vl  Vl 

Hence 

p+L  *  Vi 

1  ■*  n 

Fort-*-*,:  fa(t)  ~  +  p)t. 


fs<t)  tj 


Hence 


'ki 


Therefore 


p  ■  1/2  <k1  +  Vl>.  *  ‘/2  (Vl  -  V 

^  et,  for  g  sts  g.,  we  have  t  =  g  +  t  ,  so 
J-l  J  j-l  j 

j-l  n 

fa(t)  =Ij  ^,-1  +  '  8r>  +S  4r  '6r  *  S.  *  t.)  +  p(g.  +t  )  +  q 

r=l  J  J  r=j  r  J  J  J-l  J 
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'S«>  ;f|  W  «r-l’  +  ki‘i  *  dl  • 

where  gQ  =  0. 

Thus 

1  P 

fa(t)  =  -|>r  +  P>tj  +  <^r  "  pr  +  p)  gj-l  4 

+  i  arg-  +  f*r8r  +<l 

i  r  r  j 

fS(t>  =  V]  +Pr(8r‘gr-l)+dl 

For  t=tj  sgj: 

fa(t)  =  (p-suj)^  +^+q 

fs(t)=  kjtj+dj 

Vp-$w*4*, 

From  before: 

k^  =  P  +  f^J  P  =  1/2  (kl  +  kn+1>'  =  1/2  (k°+I  "  kl> 

j"l  n 

For  t  =  g:.,  +t. :  k  =  p  -  ?>  s 

j  1  j  j  r  j 

F°r  t=Sj-2  +  *J-1  k,-lP+'pU* 

Hence  k  -  k  ,  =  2  H 

J  J  1  J"1 

Thus,  given  k . kR+  y  gj . gp:  dj  in  the  fs(t)  form,  we 

can  determine  the  f®(t)  form  via 
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p  =  1/2  (kj  +■  kn+1) 
g  =  dj  -  1/2  (kn+1  -  kj) 

Mj  =  1/2  <kj+1  -  kj) 

Conversely,  given  the  fa(t)  form  with  Up  ....  Hn;  p.  q  and 
gj....  i6n*  we  can  determine  dp  kj,  ....  kn+p  for  the  fs(t) 


form  via 


di-q*p“j 


j-l  n 

k  *P+S^  -^s  *  + 

■’  1  J  J 


Figure  2 
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Because  we  are  dealing  with  a  piecewise  linear  and  convex 
problem,  we  have  been  able  to  reduce  the  portion  corresponding  to 
one  variable  t  to  a  separable  linear  form  on  the  t.  with  individual 

J 

interval  constraints.  Going  back  to  our  total  functional  (which  is  the 
sum  of  functions  of  one  variable  each)  and  splitting  each  variable  whose 
functional  is  nonlinear,  into  a  sum,  as  above,  we  have  thereby  reduced 
the  solution  of  the  general  convex  problem  to  that  of  a  linear  program. 
The  solution  of  the  latter  will  be  exhibited  in  the  following  section. 


4. 1  Explicit  Solution 

Because  of  the  equivalence  we  have  develop-  d  in  the  last  section 
between  functions  in  faft)  and  fs(£)  format,  we  can  evidently  always 
write  the  general  convex  piecewise  linear  separable  functional  (with 
a  finite  number  of  pieces)  as  a  linear  functional  subject  to  internal  con¬ 
straints  on  the  "incremental"  variables.  We  shall  therefore  take  our 

candidate  for  explicit  solution  in  the  reduced  (equivalent)  form: 

,a  n  ^  .(1)T  (1)  i  .(2)T  (2) 

(8. 1)  min  k  x  +  k  x 

,<1)T  (1*  j(2)T  (2)  skfe 

d  x  +  d  x 

.(1) 


(1) 


(2) 

a  s 


<  I  *  b 

A,‘”ib 


<»>  (2> 
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m  which  we  now  use  vector-matrix  notation  with 


„,1,T .  (*«« 

(i) 

(2)T  (2) 

k  s(k  j, 

n2 

(1)T  .  (1) 

X  !  (X  j, 

U). 

•  •  •  ’  xni; 

(2)T  ,  (2) 

x  =  (x  j. 

...,x(2)) 

n2 

A  =  (a. . )  is  of  full  row  rank 
ij 

Note  that  if  A  is  of  full  row  rank,  so  is  a  matrix  (J)  3  ).  Thus  the  above, 
which  we  have  written  to  separate  out  that  part  c>f  the  constraints  and 
functional  which  corresponds  to  the  non-linear  goals,  can  be  rendered 
in  the  form 

min  k  x 

(8.  2)  dTx  x  bQ 

a  <  Ax  s  b 

where  A  is  of  full  row  rank.  We  can  make  the  transformation  y  =  Ax, 
or  x  =  A  "y  +  (I-A'  K)  z,  where  A  is  a  generalized  inverse  of  A  (e.  g. ,  a 


right  inverse)  and  z  is  an  arbitrary  vector. 
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The  problem  may  now  be  rewritten 

min  kTA^y  +  k^  (I-A^A)z  with 
(9. 1)  dTA#y  +  dT  (I-A#A)  z  *  bQ 

a  i  y  *  b 


Now  let 


iTA/yj,  if  dTAj  i  0 
y.  ,  if  dTA#  =  0 


(dT  (I-A#A)k  zk  ,  if  dT(I-A#A)k  /  0 
9k=^  ,  if  (I-A#A)k  =  0 

Then  the  problem  may  be  written 


(9.2)  min  k^n  + 


«T  *T 

e.  t.  e  n  +  e  s  *  b 


where 


ifdTA#  /  0 
if  dTA^  =  0 


'l,  if  dT(I-A#A)k  t  0 
0.  if  (I-A*!k)k  =  0 
"kTA  #  /  dTAj  .  if  dTA#  t  0 

St  ,.^'  =  0 
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r*v 

""  %  rn 


A#A)r/dT(I-A#A)  if  dT(I-A#AL  ^  0 


r  ^kT(I-A#A)r  if  dT(I-A#A)r  =0 
T-# 

dAAja.  ;if  d  Aj  >  0 


'll  *j.  if^'-o 


dTA#b.,  if  dTA#  <  0 
3  3  3 


,T7#l  t#  ,T 


-# 


b  = 
3 


d  A"b  ,  if  d  A  >  0 
3  3  3 

T-# 

b  ,  if  d  A.  =0 
3  3 

,T-i  ,T -J  <  n 

d  A  a.,  if  d  A  *  0 

3  3  J 

Combining  the  components  of  n  and  s  into  a  single  vector 
T 

§  =  (Sj,  . . ^j)  in  the  order  of  decrease  of  their  functional  coefficients 

At 

kj ,  Itj, ,  and  adjoining  -® ,  +  -  lower  and  upper  bounds  for  the  ^  which  are 
Sr>  s,  the  problem  may  be  written  in  the  form 


(10.1) 


min  K  5 
with  e^§  *bQ 


a4?  sE 


where  Tfj  *■  k2  s  .  . .  s  itn,  the  ^  correspond  to  aj's  or  the  cor- 

A 

respond  to  Ex's  or  +®  the  e^  are  the  appropriate  or  For  consistency  we 
n  ^ 

must  have  5  e j  aj  s  b.  Setting  f .  =  ^  =*■  * j,  the  problem  may  be  rewritten, 
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max.  £  (-M  C, 

}  3  J  j  1  * 

(10.2)  with  £e,  j.  s  bQ  -£?.  a.  •  b 

jJJ  U  j  ^  ^  O 

“  /%,  ^ 

V)'Vki 

V 

This  is  a  linear  programming  problem.  It  therefore  has  the  following 
dual 

min.  w  b  +  £w  fj' 

0  0  j  1  1 


(10.3)  with  w  e ,  +  w.s-ki 

o  .1  j  J 


w  ,  w  so. 
°  J 


We  now  introduce 

(11)  Jq  a  {j:  -kj  *  0  and  =  o} 

J1  s  ”**j  and  ®j  = 


r- 


r  =  max  i  such  that  -k;  s  0,  e4  "  1,  Zj  e  b  s  b 

J  J  1  9  B  ( 


0  ,  / 

> 

o 


and  shall  show  that  an  optimal  solution  to  the  primal  problem  is 

•  'V/  l.' 

vv j  Jo 


(12)  5j  =  b',  j  e 


fi  -  b0  -2  e.B'  =  bi  -L 
3  °  jcJj  3  3  1  /VJo3  J 


-* 


?j  =  °<  3 6  u  J0  u  Jl 
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where  jr}  is  the  singleton  set  with  element  r.  Then 

(13)  r»-  k.  I*  =  E  (-ki)B'  +  S(-E)b'  +  (-k  )  (b'  -  £  *  > 

]  ^  J  J0  J  J i  J  .1  r  JoUJiJ 

=  (-k  )b'  ■<  Z)  f( -Tc )  -  e1(-kr))'I5:  +  terms  under  J,. 
r  o  j  J  j  1 

We  now  employ  the  regrouping  principle  and  posit  the  following  values 

for  the  dual  variables:  WQ  =  "kr,  wj  '  "kj  ”  Pj(*kr),  j  £J0  U^i  w^  -  0  otherwise 

Noting  that  these  are  all  non-negative,  we  check  that  these  w  satisfy 

the  other  dual  constraints--viz,  , 

(14. 1)  w  e  +  w  .  =  (-k  J  e  +[-k,  -  ei  ( -k_)J  =  ~k.  for  j  e  J  u  J. 

oJJ  ‘J  J  J  1  j  oi 


Thus  the  J^u  Jj  dual  constraints  are  satisfied.  Next, 

(14.2)  woe,  s  v\e  *-  -(j<r)  ej  for  jtf  JQ  J  Jj 

We  need  to  show  here  that  -(EL)  e .  a  -  k.  •  If  "k  *  0  then  e1  =  1  and  j>  r  so 

1  .1  3  ■>  J 


that  -k  «  -k..  If  -It.  <  0,  then  since  (-It  )e,  s  0  we  have  (-EL)ej  z  E;.. 
r  3  3  r  J  ~  J  3 

Thus  all  dual  conditions  are  satisfied. 

Passing  back  from  ?  to  §,  our  optimal  solution  for  (10. 1),  is  therefore 
*  ~ 


V  V  j  e  J°  UJi 


(15)  ?,  =  bn  -£  e,  b-  -L  e .  a 

r  °  JouV  3  1*  J0^  3  J 
§%aJ(  j  <5  J0  U  Jj 
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4.  2  Further  Computational  Reduction: 

For  simplicity  of  exposition  the  general  separable  convex 
piecewise  linear  functional  was  transformed  directly  into  linear  interval 
programming  format  employing  additional  variables,  one  for  each  of 
the  linear  segments.  Our  fundamental  lemma,  however,  makes  it 
possible  to  reduce  (sometimes  very  substantially)  the  number  of  seg¬ 
ments  which  need  to  be  considered  and  also  reduces  the  number 
of  variables  and  constraints  which  need  to  be  introduced.  Because 
of  this  lemma,  we  can  restrict  an  original  variable  in  the  goal  functional 
to  be  bounded  above  by  the  first  g^  at  which  the  goal  functional  minimum 
over  a  *  t  s  b  is  attained. 

The  use  of  this  technique,  together  with  other  reductions  that 
are  possible  for  special  goal  programming  functionals,  will  be  presented 
in  additional  reports  in  this  series  together  with  applications  involving 

various  utilizations  of  advance  starts.  Asa  start  toward  that  end  we  have 
supplied  a  numerical  illustration  which  illustrates  some  of  the  preceding 
developments. 


S.  Numerical  Illustration: 

We  now  conclude  the  present  paper  with  a  numerical  example 
which  illustrates  the  remarks  we  have  just  made.  For  this  purpose 
we  adapt  an  example  from  the  paper  which  initiated  this  series  of  studies 
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« 


in  manpower  planning. 

Applying  the  procedure  of  model  approximation  sketched  in 
the  next  to  last  paragraph  on  page  II -8  of  [  ]  to  a  curtailment  of  that 
example  we  have  as  our  problem 

min.  |Nj  (1)  -  301+  |n2  (1)  200 j  +  J Nj<2)  -  70|  +  |N2(2)  -  300| 

subject  to 


0  £  Njd) 

0  s 


(16) 


151^(1) 
44  £ 

143  « 


N2(l) 

+13N2<1) 


N  (2) 

15Nl(2) 


s.  <*> 


£  ® 

£  3,000 

£  a 

N2(2)  £  » 

13  Ng  (2)  £  4,  000 


Here  we  have  confined  ourselves  to  a  model  involving  only  two  types  of 
manpower  in  a  2-period  plan.  Thus,  the  variable  N^(t)  refers  to  the 
amount  of  the  i**1  type  of  manpower,  i  ■  1,  2,  scheduled  for  recruitment 
in  period  t  =  1,  2.  The  admissible  values  for  these  variables  are  in 
interval  form  except  for  the  two  budgetary  constraints  which  are 
applicable  at  the  indicated  salary  levels  of  15  and  13(thousand  dollars  /annum) 
in  each  period.  Essentially  then  this  is  all  in  the  form  of  the  expressions 
given  in  our  General  Lemma  (see  section  2,  above)  with  the  goals  of 
30,  200,  70  and  300  stipulated  in  the  absolute  value  functional  for  (16). 

By  our  Lemma  1,  the  above  example  may  be  transformed  into 
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maximize  NjU)  +  N2<1)  +  Nj(2)  +  N2(2) 
subject  to 


0  s  Nj(l) 

£  30 

0  *  N2(l) 

£  200 

15Nj(1)  +13N2(1) 

£  3,000 

44  s  Nj(2) 

*  70 

143  s 

N2(2)  ?  300 

15Nl  (2)  +13N2(2)<  4,000 


We  now  have  only  an  ordinary  linear  programming  problem  that 
involves  no  more  variables  than  the  nonlinear  problem  (16).  Further¬ 
more,  (17)  may  be  split  into  the  following  two  separate  subproblems, 
maximize  N^(l)  +  N2(l)  subject  to 


0  sNjd) 

*  30 

(18.1) 

0*  N2(l) 

£  200 

151^(1)  +13N2(1) 

*3,000 

and 

maximize  Nj(2)  +  N2<2)  subject  to 

44  *  Nj(2) 

*  70 

143  «  N2(2) 

*  300 

15Nj(2)  +13N2(2) 

£  4,000 

Without  further  ado  we  can  then  write  our  optimal  solution  for 


_ , _ _  -,.-rrril ,._.  ,, 
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(17),  and  hence  (16),  via  (18.1)  and  (18.2)  This  gives 
N.’  (1)  =  400/15,  Njd)  =  200 
Nj(2)'44  .  N*(2)=  143  4  ii|L=  257. 

Evidently  our  lemma  has  produced  a  simplification  so  that  achieving  a 
solution  to  the  nonlinear  problem  (16)  involved  practically  nothing  more 
than  merely  a  solution  by  inspection  of  two  much  smaller  linear  problems. 
This  suggests  further  possible  ways  in  which  the  ideas  contained  in  our 
lemma  may  be  extended  to  other  classes  of  problems.  Such  develop¬ 
ments  are  best  delayed,  however,  until  they  can  be  treated  in  their 


own  right  via  other  papers  in  this  series. 


-  (j; 
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